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Abstract 


We analyze the approximately 150 reported gamma burst 
spectra of the Konus catalogue (Mazetp e^ al.. 1981) using the 
thermal synchrotron model. An overwhelming majority of these 
spectra can be satsifactorily fitted by theoretical thermal 
synchrotron spectra of mildly relativistic electrons in strong 
magnetic fields, making the strong-field neutron star picture at 
least self-consistent. Valuable additional information is also 
extracted from various spectral features contained in many of the 
events. 


1 



I • INTRODUCTION 


In a recent article we pointed out that many of the observed 
gamma-ray burst spectra can be fitted by optically thin thermal 
synchrotron (or high -harmonic cyclotron, hereafter abbreviated 
TS) emission of mildly relativistic (kTjj ~ mc^) electrons (Liang 
1982). In most cases the TS fits are at least as satisfactory as 
the conventionally adopted thermal bremstralung (TB) or free-free 
spectra. The success of the TS fits corroborates the picture of 
strongly magnetized neutron stars as gamma-burst sources. Since 
then we have performed a detailed and systematic analysis of the 
entire Konus Catalogue (Mazets e^ aX. 1981) containing spectral 
data of ~ 120 separate events and 150 burst spectra (20-1500 keV) 
obtained by the Konus experiments from 1978 through 1980. This 
article summarizes the major findings. Soms of these results 
have far-reaching implications for future model building of the 
gamma-burst sources. 

Despite the controversy surrounding the alleged emission and 
absorption features reported by Mazets ^ Jt a.X» (1981) in many of 
the events (due to uncertainties introduced by the deconvolution 
procedure [Fenimore et. aX* I 982cl), we believe that the overall 
continuum shapes of the spectral data remain trustworthy, since 
similar continua are also observed by other groups (e.g., Cline 
e t a 1 . 1975, Gilman ^ 1980). But it remains interesting for 

the experimentalist to try to deconvolve Mazets' data with an 
assumed TS instead of TB spectrum to see if there are significant 
d if f erences . 
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IT 


SPECTRAL FITS 


Using the analytic formula derived by Petrosian (HUl) which 
is valid in the regime^' >> ~ eB/2?rroc is the Lamor 
frequency and T ^ RTg/mc^; see also Trubnikov 1961, Lamb and 
Masters 1979), we attempt to fit all 150 spectra of Mazets' 
catalogue wicn the TS model. Since chi-square analyses of the 
entire collection is prohibitively labor consuming, we simply 
used variable-endpoint fitting and eyeball inspections for the 
goodness of the fit. Overall, we are convinced that over 80 
percent of the fits should be considered satisfactory. The 
remaining spectra either have very poor statistics or have shapes 
clearly distinct from the monotonic convex shapes of TS (or TB) 
spectra. These anofflaloijs spectra, however, may still be produced 
by cyclo-synchrotron emission of nont h erm al electrons, or they 
could simply be the composite of time-varying spectra due to 
changing electron distributions within the ~ 4 sec integration 
period. In fact, considering the ultra-short cooling time of the 

electrons (tcooling ‘^Couloumb relax* independent of whether 
the cooling is due to TS, TB or inverse Compton (IC) as proposed 
by Fenimore ejt ^JL . [1982a]), it is just remarkable that most 
sources can maintain close to thermal distributions at all. 
Collective plasma processes must be extremely efficient in 
thcrmalizing the electrons, a clue to be noted by all model 
builders of gamma-ray bursts. 

Even among the "good” fits, there are systematic trends of 
small deviations from the TS continuum shape. Some of the more 
significant ones are listed here. (1) Many of the spectra 
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deviate from the TS fit at the high~6nergy tail (moatly dropping 
below). This may be due to systematic error (e.g,, detector 
deficiency near detection limit). If it is intrinsic, it would 
suggest that the h igh-ene r g y- 1 a i 1 electrons often fail to stay 
Maxwellian, an expected phenomenon since they have the fastest 
cooling and slowest relaxation rates. A few of the spectra with 
high-energy excess may also be due to s elf-Compton izat ion (Liang 
1981). (2) Some spectra have apparent excess emission in mid- 
energies (100-500 keV), If they are real, it could be due to 
broadened emission lines (e.g., redshifted 511 keV line or 
cyclotron line with B >, 10^^ G) or Comptonization. (3) Some 
spectra show low-energy cutoff. Those would most likely be due 
to self-absorption. (4) A number of spectra show conspicuous 
excess emission at the low-energy end. In the TS picture, it 
would be most natural to attribute them to low-energy harmonic 
emissions (cyclotron lines). Several spectra, in fact, show 
double peaks characteristic of the first and second harmonics 
(Sec, V). The undisputed detection of the first few harmonics 
would be the strongest evidence for the TS model. (5) Seme 
spectra have "absorption" dips in the 30-70 keV rangCi but the 
continuum continues to rise below the frequency of the dip. The 
reality of these absorption features has been challenged by some 
(Fenimore ^ aJL. 1 982c, 1982b). Even f they are real, their 
interpretation is by no means obvious, since the rise of the 
continuum below the dip energy might suggest that the field 
strength responsible for the TS continuum is weaker than the 
absorption region. One way out is that since the TS continuum 
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emissiun region is much hotter than the line absorption region 
(as evidenced by the narrowness of some of the lines), the first 
harmonic emission peak would be Doppler shifted down (v £ » 

+ Cy T) c^ ~ 2-3) to undetectable energies, while the absorption 
feature would stay near the unshifted value 

Figures 1-6 provide samples of TS fits to the data and 
representatives of the above five categories of deviants. 
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III. DISTRIBUTION OF ^ 

The spectral fitting uniquely determines the ’’critical 

frequency” « Vj^ T^ <sin 0>, where 0 is. the angle between 

magnetic field lines and the line of sight. In Figure 7a we plot 

the histogram for the distribution of Vg for all the ~ 150 

spectra. A more meaningful histogram is the distribution of 

in terras of events (■* 120), where we pick only the hiRhest v ^ for 

each event with reported spectral evolution. It is clear that 

the distribution peaks between 2-5 keV and cuts off sharply at 

~ 12 keV. It we assume a universal field strength of B<sin 0> “ 

2 X 10^^ G = 23 keV, clearly not true in some cases), then 

tne temperatui’e distribution peaks between .2-. 7 mc^ (top scales 

2 

of Fig. 6). It is significant that few events have T >. me 
(unless B<sin 0> was << 10^^ gauss). 


IV. TIME EVOLUTIONS 


Mazets' catalogue included data on the spectral time- 
evolution for about 15 of the events. In Figure 8 we plot the 
time development of for 14 of these events (the other one has 
a late-time spectrum deviating so badly from TS that no 
meaningful could be determined). They show a distinct trend 
of initial exponential decay (with decay constant “ sec) but 
leveling or oven bouncing back after a couple of e-foldings. 
This should be a useful clue to modeling the energy source. Two 
ot the events also show evidence of eorJy-time self-absorption 
and late-time harmonic emissions. We will discuss the important 
implications of such features for these two events in a separate 
article (Liang 1983). Table 1 tabulates the crudely estimated 
decay constants for v ^ of all events. 


V, 


SPECTRA WITH HARMONIC FEATURES 


Taken at face value, ” 24 of Mazets' spectra show 
discernible low-energy features (in excess of the TS continuun), 
inter pretable as first harmonic em^\sion. Those spectra allow us 
to determine B and T separately (Table XI). It is interesting to 
note that, overall, the ones with observable harmonics tend to 
have higher v^, than the overall norm. Eight of these events also 
exhibit possible second harmonic peaks (at ZVjl). The ratio of 
the first to second peak, as well as upper bounds to the line 
width are also consistent with the theoretically predicted 
values, except for the event of 5-24-79 (see Table III). In 
Figure 9 we plot the histograms for the T and B distributions. 
It is remarkable how narrow the distributions are. We believe 
the existence of harmonic features is the strongest evidence in 


favor of the TS model. 


VI. SPECTRA EXHIBITING SELF-ABSORPTION 


Approximately a dozen ovonta contain apectra. \rith low-energy 
cutoff intorpretab le aa a e I f -obs or p t ion. For these spectra wo 
can determine n^h (electron column dcnsity)/TK 2 (T”^) (K 2 Bessel 
function, cf. Eq. 13, Liang 1982). If, moreover, we assume that 
vl lies at or below the peak, then lower bounds to T, n^h, L^y^^/A 
(synchrotron luminosity per unit emission area, the coefficients 
in Eqs. (11), (14), etc., of Liang 1982 should be lowered by a 
factor of ^ 3) and d (source d i s t anc e) / A^ ^ ^ can be estimated. 
Table IV lists these values. Despite the uncertainties in the 
distance estimates, some of these events are evidently extra- 
galactic (unless A is tiny), s^iowing that the 3-5-79 event is not 
unique (Raraaty a,JL. 19 81, Liang 1981), We must emphasize, 
however, that these few sources likely represent the most distant 
ones . since all others do not show se If -absorpt ion down t 15-30 
keV. Hence the TS interpretation is con s istent with othe.; evi- 
dence showing that the vast majority of gamma-burst sources are 

galactic. Also, since TK 2 likely lies in the range 0.5-0.05, the 

21 2 2 — 2 

electron column densities are uniformly around 10 -10 cm 

These are remarkably optically thin sources. 
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VIX. COMPARISON BETt^' E^N TSj. TB, MJ) ^ SPECTRAL FITS 

For a spectral fit that is geared towards the mid-energy 
data points (50-500 keV), the three theoretical spectra, namely, 
TS , TB and IC should in principle be distinguishable at the low- 
and high-energy ends. Unfortunately, most spectra do not have 
enough data points at the extremes to distinguish clearly among 
the three models. Moreover, low-energy absorption and high- 
energy non-thermai emissions may also mask the underlying thermal 
spectra. Figure 8 represents one of those "clean" thermal 
spectro with a monotonic curvature extending over the widest 
observed energy range (20-1500 keV). Plotted on it are three 
best-fit theoretical spectra corresponding to TS, TB and IC, It 
is clear that TS and TB are almost indistinguishable lexcept at 
tne highest energy, where the data are poor. On the other hand, 
IC is not so good a fit overall despite its three free parameters 
(compared with two for the others). The high n^h required by the 
IC fit and the high n^h of the TB fit would make them 
incompatible with the strong-field neutron star picture (Table 


vm. 


SUMMARY 


Xt ie cleor that a huge amount o£ uaeful information it 
hidden in the data publiahed by Mazets (1981). We are 
convinced that an overwhelming majority of these spectra can be 
satisfactorily interpreted as thermal synchrotron radiation with 
^ X 10^® G) T^) lying between 1-12 keV, although this 
interpretation is not unique. This is basically consistent with 
the picture of hot (kT ** mc^) thermal emission from strong field 
(B “ 10^^ G) neutron star surfaces. Various spectral features 
provide additional valuable information about the source, once we 
adopt the TS interpretation. In particular, several sources 
showing 8 e 1 f - ab s 0 r p t io n can be estimated to lie as far as the 
3-5-79 event (e.g,, the Magellanic clouds), and sources with 
detectable first harmonics have vary narrow distributions in 

i ; ur e and magnetic field strength. 

Clearly, we have only scratched the tip of the iceberg. 
Future analyses taking into account, e.g., the correlation 
between spectra and intensity variations, annihilation-line 
features, etc., promise to provide even more information. It is 
to be hoped that such information will ultimately suffice to 
reveal the true nature of the majority of the events, which 
likely have a rather universal origin. Ruderman's wishes (1975) 
may yet be fulfilled after all. 
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Table I . Initial Decay Time of ^ 


Event 

a t/A In V )/ sec 

9-14-78 

6.9 

9-18-78 

-9.6 -> +9. 

11-04-78 

2.3 

11-19-78 

5.8 

11-24-78 

2.4 

3-05-79 

6.0 

3-13-79 

5.7 

3-25-79 

13.2 

3-29-79 

4.1 

4-02-79B 

5.9 

6-04-79 

3.7 

10-14-79 

2.3 

11-09-79 

3.2 

11-11-79 

9.3 
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Table II. Parameters of Spectra with First Harmonic Emission 


Spec.trum^^^ 

V ° ( k e V ) ^ ^ 

Vg (keV) 

T(kTg/mc^) 

J(10^^ gauss) 

9-16-78 

40 

11 

.36 

7.32 

9-18-78^ 

20 

1.8 

.23 

2.93 

10-12-78B 

43 

8.8 

.32 

7.27 

11-09-78 

35 

4.1 

.26 

5.23 

1-19-79 

35 

3.4 

.24 

5.09 

2-13-79 

28 

5.3 

.32 

4.46 

2-15-79 

28 

9.5 

.39 

5.38 

3-13-79^ 

32 

10 

.38 

5.97 

3-13-79^ 

32 

0.6 

.12 

3.59 

3-29-79^ 

18 

10 

.49 

3.59 

3-29-79® 

18 

0.032 

.04 

1.72 

4-02-79B® 

32 

13 

.43 

6.32 

4-24-79A 

18 

2.7 

.29 

2.77 

5-24-79 

18 

5.0 

.37 

3.15 

5-25-79 

50 

3.8 

.22 

6.77 

6-12-79 

32 

3 .4 

.25 

4.69 

6-28-79A 

32 

7.4 

.34 

5.52 

7-05-79 

32 

1.6 

.18 

4.26 

7-12-79 

23 

1.9 

.22 

3.38 

10-07-79 

18 

3.3 

.31 

2.96 

11-09-79^ 

30 

1.8 

.19 

4.30 

11-11-79® 

30 

2.5 

.22 

4.29 

1-16-80 

30 

3.0 

.24 

4 .45 

1-26-80 

40 

1 .4 

.16 

5.10 
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[a] Superacripts above event date denote time interval of spec- 
trum for events with reported sp.ectraX evolution (see 
Mazets e,jt a_l.. 1981). 

lb] denotes the estimated peak of the obs erved emission 
feature. Since in most cases there is no data with energy 
below v£, the actual first harmonic peak may lie slightly 
below the listed value. In that case the actual B would be 
slightly below, and the actual X would be slightly above the 
listed values. Because of these uncertainties, the apparent 
evolution of B-field for the 3-l'3-79 and 3-29-79 events should not 
be taken at face value. 

l c] T - '^c ^ r ( T ) where the relativistic Doppler factor 

['(TJ rr 1 + c T wihh c ~ 'n 

ranges (see numerical results in Lamb 1982). No allowance 
for gravitational redshifting (which could be as much as 20- 
30 percent) has been included in these T and B estimates. 
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Table 111. spjg>ctj^ Ikii M± Iftcoad 


Hnyroonic Em ia a Ion Feature* 


Spectrum 

1st Harmonic 
Location 
\>X (keV) 

2nd Harmonic 
Location 
\>2 (keV) 

^ 2 
Flux Ratio 

Observed 

TS 

Prediction 

9-16-78 


40 

80 

1.25 

? ,32 

10-12-78B 


43 

85 

1.83 

1.57 

11-09-78 


35 

70 

2.0 

1.92 

4-02-7 9 


32 

65 

1.32 

1.16 

4-24-7 9A 


18 

36 

1.85 

1.72 

5-24-79 


18 

36 

4.0 

1.35 

6-28-7 9A 


32 

65 

1.36 

1.47 

7-05-79^ 


32 

65 

2.5 

2.78 

Note: 

In some 

cases the 

omission features 

can in principle 

be viewed 


cyclotron absorption features of lower field strength. But the excellent 
agreement between columns 4 and 5 in almost all cases suggests that the 
harmonic emission interpretation is likely correct. 
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Table IV. Spectra vith Self"AbBorption 


Spectrum 

t 

(keV) 

^c 

(keV) 

ngh/TK2(T‘"l) 

(10^^/cm^) 


Assumingv j 

li'^ab 

• 

1 

Te.ClO-'^) > 

L > 

syn — 

d(kpc)* > 

A(km^) 

A^^^Ckm) 

10-06-78A 

30 

9.9 

i.6 

.39 

2.2 

4.5x10^® 

30 

11-04-78A 

25 

2.1 

5.2 

.22 

.34 

3.1x10^^ 

10 

U-07-78 

42 

6.0 

4.8 

.29 

1.7 

3.3x10^® 

37 

11-11-78 

28 

1.9 

7.7 

.20 

.24 

3.2x10^^ 

11 

11-19-78^ 

40 

110 

0.6 

.87 

.20 

4.9x10^^ 

38 

11-19-78^ 

40 

6.0 

4.3 

.29 

1.5 

3.0x10^® 

10 

11-21-78A 

26 

4.8 

2.3 

.31 

1.2 

1.1x10^® 

19 

2-14-79 

A 

JU 

4.1 


no 
• <c.o 

1 A 
Oi • V 

1 i«in40 

A • A4WA W 

1 *1 

A*/ 

4-02-79B^ 

60 

50 

1.7 

.56 

9.7 

1.7x10^2 

49 

4-06-79 

45 

38 

1.2 

.56 

6.9 

7.5x10^^ 

50 

6-13-79 

50 

1.9 

32.2 

.17 

.33 

9.3x10^^ 

72 

10-14-79^ 

40 

5.9 

4.4 

.28 

1.1 

2.8x10^° 

15 

11-11-79^ 

50 

5.9 

6.9 

.24 

0.6 

3.8x10^*^ 

30 


^The self-absorption frequency is approximately where the optically thin 
spectrum would intersect the Rayleigh- Jeans part of the black-body spectrum. 


*Due to anisotropy of the radiation pattern and unknown B-field orientation 
ana other factors, these distance limits must be considered as uncertain by at 
least a factor of 2. Moreover, no gravitational redshifting has been included 
here, which would lower d by an additional (1 + 
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Table V. Comoar Ison Between Parameters Determined by 

TS . TB and IC Pita to the Spectrum of 9-18-78 ^ 



TS (Vj, - 6.2 keV) 

TB 

IC 

T 

(assume v£ < 25 keV) 
> 180 keV 

317 keV 

200 keV 

^ es 

(assume < 25 keV) 

< J . 5 X 10“^ 


"I 

n^h 


(assume < 25 keV) 

1 10^^ 
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FIGURE CAPTIONS 


Figure I 


Figure 2 


Figure 3 


Representative gamma-burst spectra typical of the 
"pertect fits" with the thermal synchrotron model. 
Values of corresponding to the tneory cur 'fee ore 
listed in upper righthand corners. 

Representative spectra with high-energy tail 
deficiencies compared to the theoretical TS curves. 

Representative spectra with emission excesses above the 
TS curves in the 300-800 keV range. (a) and (b) are 
typical of "broad" features, while (d) is typical of 
the narrow "line" at ** 400 keV, interpretable as the 
redshifted 511 keV line. (c) has the unique property 
that the feoture was very broad at early times (spec- 
trum 1, observation by U.S. groups, may have resolved 
it into two separate lines at 400 and 800 keV), but 
evolved into a narrrow line at ~ 600 keV at later times 
(spectrum 3), which, if real, can correspond only to a 
blueshifted 511 keV line. The feature is not evident 
at time 2. Note that spectra 1 and 3, which have the 
emission feature, also show evidence of self-absorption 
at the low-energy end, probably corresponding to higher 
total electron (pair) column density due to the abun- 
dance of pairs. 
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Figure 


F igure 


Figure 


Figure 


Figure 


Figure 


Figure 


4 Representative spectra exhibiting aelf-alsorption. 
Dashed curves are poasibLc Rayleigh -Jeans limit s 
corresponding co the appropriate electron temperature. 

5 Representative spectra exhibiting first (or even the 
second) harmonic omissions at energies below the TS 
continuum. Dashed curves suggest possible profiles of 
the harmonics and are not detailed model fitting. 

6 Representative spectra with the controversial cyclotron 
absorption dips at 40-65 keV. Note that the TS fits 
continue to be excellent below the dips. 

7 Histogram of distribution according to the total 

number of published spectra (a) and events (b). 

8 Time evolution of ^ for 14 events with reported 
spectral evolution. 

9 Histograms of T and B distributions for the 22 events 
with possible first harmonic emission features. 

10 Comparisons between TS , TB and IC fits to the continuum 
spectrum of 9-18-78.^ The narrow feature at 400 keV is 
removed from the fitting procedures. 
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